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Abstract 



(N 
Ph. 

Ph' We present a supersymmetric T13 model for the tri-bimaximal neutrino mixing, and 

r~| ' the complete flavor group is T13 x Z4 x Z2. At leading order, the residual symmetry 

of the charged lepton sector is Z3, and the T13 symmetry is broken completely in the 
CN ! neutrino sector. The charged lepton mass hierarchies are determined by the sponta- 

I neous breaking of the flavor symmetry, both the type I see-saw mechanism and the 

I Weinberg operator contribute to generating the light neutrino masses. Tri-bimaximal 

mixing is exact at leading order while subleading contributions introduce corrections 
of order to the three lepton mixing angles. The vacuum alignment and subleading 
■ corrections are studied in detail, a moderate hierarchy of order Ac between the vacuum 

expectation values of the flavon fields in the charged lepton and neutrino sectors can 
be accommodated. 

X 



^Email: dinggj@ustc.edu. en 



1 Introduction 



The presence of two large and one small mixing angles in the lepton sector [TH3], 

0.27 < sin^ 612 < 0.37, 0.39 < sin^ 623 < 0.64, sin^ < 0.040 (0.044) at 3a (1) 

suggests that the observed neutrino mixing matrix is remarkably compatible with the so 
called tri-bimaximal (TB) structure [1] within measurement errors. The simple form of the 
TB mixing matrix implies an underlying family symmetry between the three generations 
of leptons. It has been realized that the TB mixing can naturally arise as the result of a 
particular vacuum alignment of scalars that break spontaneously certain flavor symmetries. 
In the past years, much effort has been devoted to produce TB mixing based on some family 
symmetry. It has been shown that TB mixing can be understood with the help of discrete 
flavor symmetries, such as A4 [514T0]. T7 [11] , T' [12], S4 [T3l[T^ and A(27) [I5], or continuous 
flavor symmetry 50(3) [IS] and SU{3) [T7]. Discrete non-abehan groups appear to be partic- 
ularly suitable to reproduce the TB mixing pattern, some higher order discrete groups such 
as A5 [18], A(54) [19], S(81) [20] and PSL2{7) [21] are also considered for neutrino mixing 
besides the above mentioned simple groups, the extension of the discrete flavor symmetry to 
the quark sector and grand unification theory (GUT) have been investigated as well [TJIH], 
please see Refs. [221123] for a review. In this work, we shall study another 39 element simple 
discrete group T13 in flavor model building, which has gained much less attention. 

Recently 76 discrete groups with 3-dimensional representation were scanned, it is sug- 
gested that Ti3 is the group with the largest fraction of TB mixing models [21]. But the 
authors set all the couplings to be equal to 1, the vacuum expectation values are chosen 
to be or 1, and the vacuum alignment is not considered dynamically in Ref. [21]. It is 
very interesting to investigate the possible consistent realizations of TB mixing based on 
Ti3 group from this point of view. As far as we know, the T13 group as a discrete flavor 
symmetry has not been discussed extensively. We note that a T13 flavor model was put 
forward in Ref. [2S], and its implication in the indirect detection of dark matter was studied. 
However, the motivation is not to produce the TB mixing B We have tried many possible 
assignments for the involved fields, we find the TB mixing can be produced exactly at leading 
order (LO) in some scenarios, but meanwhile we face the difficulties that the first and third 
light neutrino are degenerate or the corresponding vacuum alignment is very difficult to be 
realized or some other problems. In particularly, the realizations of TB mixing based on T13 
symmetry are drastically constrained after taking into account the vacuum alignment issue. 
After lots of trial and error, we construct a T13 flavor model described in this work, where 
TB mixing is obtained exactly at LO. It is well-known that discrete group or continuous 
one like U{1) are usually introduced to eliminate unwanted couplings, to ensure the need 
vacuum alignment and to reproduce the observed charged charged lepton mass hierarchies. 
In the present work, the auxiliary symmetry Z4 x Z2 is introduced for this purpose. It is 
notable that the charged lepton mass hierarchies are determined by the T13 x Z4X Z2 flavor 
symmetry itself without invoking a Froggatt-Nielsen U{1) symmetry. 

This paper is organized as follows. In section 2, we discuss the relevant features of T13 
group. In section 3, the structure of the model is described, the LO results for neutrino as well 

^The vacuum alignment and the next leading order correction are not discussed in Ref. H5], a set of 
numerical values are chosen by hand for the model parameters so that the resulting lepton masses and flavor 
mixing arc consistent with experimental data. 
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as charged lepton mass matrices are presented. In section 4, we show how to get in a natural 
way the required vacuum ahgnment used throughout the paper. In section 5, we present 
the study on the corrections introduced by the higher order terms, which is responsible 
for the deviation from TB mixing. Finally section 6 is devoted to our conclusion. We 
give the explicit representation matrices and the Clebsch-Gordan coefficients of T13 group in 
Appendix A. The analysis of the subleading corrections to the vacuum alignment is presented 
in Appendix B. 



2 The discrete group T] 



13 

The discrete group T13 is a subgroup of SU{3), and it is smallest discrete group with two 
complex irreducible three-dimensional representations. T13 is isomorphic to Z13 x |26|l27]. 
consequently it has 39 group elements. T13 can be generated by two elements S and T 
obeying the relations 

^13 = = ST = TS^ (2) 

The 39 elements of the group belong to 7 conjugate classes and are generated from S and T 
as follows, 



C2 
C3 

C7 



T, TS, TS^, TS\ TS\ TS^ TS'^, TS^\ TS^\ TS^^ 

rp2 j^^g T'^S'^ T'^S^ T^S"^ T'^S^ T'^S^ T^S"^ T^S^ T^S^ 2^2^10 2^2^11 rp2^i2 

s, 5^ 

^10 ^12 



7 c8 cll 



S\S\S' 



(3) 



The Ti3 group has 7 inequivalent irreducible representations li, I2, I3, 3i, 3i, 82 and 82. 
It is easy to see that the one-dimensional representations are given by 



li 
I2 
I3 



S 
S 

s 



1, T 
1, T 
1, T 



1 

u 



(4) 



where to = e^*'^/^. The three-dimensional representations are given by 

3i : S = 



S 




T 




(5) 



where p = e^*'^/^^, the 3i and 82 representations can be obtained by performing the complex 
conjugation of 3i and 82 respectively. We can straightforwardly calculate the character 
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Table 1: Character table of the T13 group, where ^1 ~ p + + p^, ^2 ~ P^ + P^ + P^, P = g2i7r/i3 g^^^^j 
denotes the number of the elements contamed in the class d, and hd is the order of the 

elements of d. 



table of Ti3, which is shown in Table [TJ Then the multiplication rules between various 
representations follow immediately, 

ll® R = R^ll = R, I2®l2 = l3, l2®l3 = ll, I3®l3 = l2, 

li®3i = 3i, li®3i = 3i, li®32 = 32, li®32 = 32, 
3i O 3i = 3i5 © 3iA © 32, 3i ® 3i = ll © I2 © I3 © 32 © 32, 
3i O 32 = 3i © 32 © 32, 3i © 32 = 3i © 3i © 32, 
3i © 3i = 3i5 © 3iA © 32, 3i © 32 = 3i © 3i © 32, 

31 © 32 = 3i © 32 © 32, 32 © 32 = 3i © 325 © 32vi, 

32 © 32 = ll © I2 © I3 © 3i © 3i, 32 © 32 = 3i © 325 © ^2A (6) 

where the indices i = 2, 3, R indicates any T13 irreducible representation, and the subscript 
S and A denote symmetric and anti-symmetric products respectively. The explicit represen- 
tation matrices of the group elements for the three dimensional irreducible representations 
are listed in Appendix A. From these representation matrices, one can directly calculate the 
Clebsch-Gordan coefficients for the decomposition of the product representations, which are 
given in Appendix A as well. 



3 The structure of the model 

The model is supersymmetric and based on the discrete symmetry T13 x Z4 x Z2. Super- 
symmetry (SUSY) is introduced in order to simplify the discussion of the vacuum alignment. 
All the fields of the model, together with their transformation properties under the fiavor 
group, are listed in Table [2l We assign the 3 generation of left-handed lepton doublets £ to 
be the 3i representation, while the right-handed charged lepton e*^, fi'^ and t'^ transform as 
ll, I2 and I3 respectively. It is notable that the three right-handed neutrinos u^, 1/3 and 
are assigned as li, I2 and I3 as well, they transform in the same way as the right-handed 
charged lepton fields. This is an interesting feature of the model. We note that in popular 



3 



Fields 


i 

















y 

A, 






T 


r; 
1 


A, 


r 


^0 


„o 

1 


^0 




3i 


li 


I2 


I3 


li 






ll 


3i 


3i 


3i 


32 


32 


I2 


I3 


32 


ll 


z. 


1 


i 


-1 


-i 


1 


1 


1 


1 


i 


i 


1 


1 


-1 


-1 


-1 


1 


-i 


Z2 


1 


1 


1 


1 


-1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


1 


1 


1 


-1 


U{1)r 


1 


1 


1 


1 


1 


1 


1 

















2 


2 


2 


2 


2 



Table 2: The transformation properties of the matter fields, the electroweak Higgs doublets, the flavon fields 
and the driving fields under the flavor symmetry T13 x Z4 x Z2. 



A4 and 5*4 models, the right-handed neutrinos are frequently treated to be a triplet 
Lepton masses and mixing arise from the spontaneous breaking of the flavor symmetry by 
means of the flavon fields, they are neutral under the standard model gauge group and are 
divided into two sets = {x^^} and $y = {0, 77}. We note that all the flavon fields are 
triplets under T13 in this work, is responsible for the charged lepton masses and $j/ for 
the neutrino masses at LO. In the following, we shall discuss the LO predictions for fermion 
masses and flavor mixings. For the time being we assume that the scalar components of the 
flavon fields acquire vacuum expectation values (VEV) according to the following scheme. 



{x)=vA 1 , (0 





{<P)=vA 1 , {V)=\ V, \ (7) 



In section 4 we shall show that the above alignment is indeed naturally realized at LO from 
the most general superpotential allowed by the symmetry of the model. 

3.1 Charged leptons 

The charged lepton masses are described by the following superpotential 

5 



i=l 

where 

o = {{xixxhjs^, ((xx)3is03is' ((xx)3is03iA' iixxh^Os,, {x{^0s2hi} (9) 

We note that the subscripts li, I2, I3, 3i etc denote the T13 contractions. In the above 
superpotential wg, for each charged lepton, only the lowest order operators in the expansion 
in powers of 1/A are displayed explicitly. Dots stand for higher dimensional operators which 
will be discussed later. It is remarkable that the Z4 symmetry imposes different powers of x 
and ^ for the electron, muon and tauon terms, i.e., only the tau mass is generated at LO, the 
muon and the electron masses are generated by high order contributions. After electroweak 
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and flavor symmetry breaking, we have 



A3 



64 



A3 



^£5 



A3 



/ c 2 

VdC^ie + /i + r) + 2y^^VdH''{e + uj h + ujt) 



e + /i + r) + y^j^Vdli 



1) 



(10) 



where Vd = {hd), the parameters ye and are parameterized as ye = yei + (4?/e2 +^64" 
ye^v'^/v'^ and = 2?/|^. As a result, the charged lepton mass matrix has the form 



me 



1 








v 




UJ 



Vd 



Obviously the charged lepton mass matrix is diagonalized by performing the transformation 
i — )■ Uei, where Ue is 

1 
1 
1 




1 1 



UJ 



(12) 



and the charged lepton masses are given by 

,3 



me 



yjf.vd 



m,. 



V3 



m-r 



V3 



yr^Vd 



(13) 



We see that the charged lepton mass hierarchies are generate by the spontaneous breaking 
of the flavor symmetry. To estimate the order of magnitudes of and v^, we can use the 
experimental data on the ratios of charged lepton masses. Assuming that the coefficients ye, 
y^ and yr are of C(l), we have 



me 



v: 



— ~ -4 ^ 0.0003 

m^ h.^ 



m 



m-r 



0.06 



These relations are satisfied for 



:^,^) ~ (0.017, ±0.032) 
' A A 



(14) 



(15) 



we see that the amplitudes of both v-^/A and v^/ A are roughly of the same order about A^, 
where Ac is the Cabibbo angle. It is interesting to investigate the fiavor symmetry breaking 
pattern in the charged lepton sector, it is induced by the VEVs of x ^i-nd ^ at LO. Given 
the explicit representation matrices listed in Appendix A, it is obvious that the VEVs of 
X and ^ are invariant under the action of T and T^. Furthermore, we can check that the 
hermitian matrix mjm^ is invariant under both T and as well. Therefore we conclude 
that the T13 fiavor symmetry is broken down to the Z3 subgroup generated by the element 
T in the charged lepton sector at LO. 
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3.2 Neutrinos 



The superpotential for the neutrino sector can be written as 

+ ^((«M3.(0r/)3ji, + ... (16) 

where Mi and M2 are constants with dimension of mass, they are naturally of the same 
order as the cutoff scale A, and the factor of | is a normalization factor for convenience. 
We note that w^^ denotes the lagrangian of the type I see-saw mechanism, and wl^^ is 
the collection of higher dimensional Weinberg operators. Taking into account the vacuum 
alignment shown in Eq.(j7]), we can read the Dirac and Majorana mass matrices immediately 
from w^^ as follows 

/ Vui -Vui \ ^ 
mD = I uj'^yu2 -^^yu2 ^Vu, rriM = | M2 | (17) 
V uy^-i -u'^yus I 

where Vu is the vacuum expectation value of the Higgs field hu- It is remarkable that the 
eigenvalues of the Majorana mass matrix mu are Mi, M2 and — M2, two of the right handed 
neutrinos are degenerate at LO. This is a distinguished feature of our model from the previous 
flavor models in which the right-handed neutrinos are usually treated to form a triplet. It 
is very interesting to discuss the assignment of right-handed neutrinos as singlets and the 
corresponding phenomenological implications in flavor models based on ^4, A(27), 6*4 and 
so on. The light neutrino mass matrix from see-saw mechanism is given by the well-known 
see- saw formula 



ml^ = -m^mJ^mD = \ -a - 2b a-b \ ^ (18) 
where 

a = yl,vl/{AMi), b = y^^yuzvH {KM^) (19) 
The superpotential ^ leads to the following effective light neutrino mass matrix 






ml" = \ s t \ ^ (20) 
where 

r = -2x^4V^v^/A'^ 
s = 2x^^vl/K^ 

t = -Ax,,vl/A' + 2x,2uJ/A' (21) 
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Therefore in the flavor basis where the charged lepton mass matrix is diagonal, the light 
neutrino mass matrices read 




-26 b b 
Tnt^ =16 a —a - 
b —a — b a 

r + 2s + 2t r — s — t r — s — t 
nfjf=[ r-s-t r-s + 2t r + 2s - 1 (22) 
r — s — t r + 2s — t r — s + 2t 



Both the light neutrino mass matrices m^^ and ml/^ are 2^3 invariant, and they satisfy 

• i. I SS(eff)^ , / S5(e//)N / SS(eff)^ , / SS(eff)^ m, r 

the magic symmetry [mu )ii + [mu jis = {'mu )22 + {i^u )32- iheretore 
they are exactly diagonalized by the TB mixing matrix 

UTBmfUTB = diag(-36, 0, 2a + 6)^ 

UTBm'JfUTB = diag(s + t, r, s + t)^ (23) 

We note that the contribution m^^ from the see-saw mechanism is of the same order as ml-^-^ 
coming from the Weinberg operators, consequently both contributions should be included. 
The light neutrino mass matrix is the sum of m^^ and ml-^-^ 

= ml^ + nfj^ (24) 

Obviously my is still diagonalized by the TB mixing matrix, and the light neutrino masses 
are given by 

m^i = {s + t-3b)-^ 

m,3 = i-s + t + 2a + b)^ (25) 



where Utb is the well-known TB mixing matrix 



U' 



TB 







/ /2 J__ 

3 

111 
\/3 



(26) 



We note that the contributions proportional to a and b can be absorbed into s and t by redef- 
inition s — )■ s — a — 26 and t ^ t + a — b, therefore the light neutrino masses depend on three 
unrelated complex parameters. There are more freedoms to tune the mass differences and 
then satisfy the constraints associated to neutrino oscillation, the neutrino mass spectrum 
can be normal hierarchy or inverted hierarchy. In contrast with some "constrained" flavor 
models, no neutrino mass sum-rules [2H] can be found in this model. We could certainly 
remove the right-handed neutrinos from our model, then the neutrino masses are described 
by the Weinberg operators wf,^^, the above conclusions remain invariant. However, if we only 
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concentrate on the see-saw realization w^^, the second neutrino would be massless although 
the lepton mixing is of TB form, this scenario is ruled out by the experimental observations. 

It is notable that the VEVs of </> and t] are always changed under the action of any T13 
group element except unit element, consequently the flavor symmetry T13 is broken down to 
nothing in the neutrino sector at LO. Reminding that ones usually break the flavor symmetry 
into the low energy neutrino symmetry group Klein four [29H3T] or Z2 [32l|33] to guarantee 
TB mixing for neutrinos, it is really amazing we can still obtain TB mixing even if the flavor 
symmetry is broken completely in the neutrino sector at LO. 

In short summary, at the LO the T13 flavor symmetry is broken down to subgroup 
and nothing in the charged lepton and neutrino sectors respectively. This breaking chain 
lets us to find the TB scheme at LO as the lepton mixing matrix. However, the mixing 
angles generally deviate from the TB values after the corrections of the higher order terms 
are included. It is remarkable that this symmetry breaking pattern of our model has not 
been studied, as far as we know. It is attractive to investigate whether we can still reproduce 
TB mixing in models with ^4 or 5*4 symmetry, if the flavor symmetry is broken completely 
in the neutrino sector at LO. 



4 Vacuum alignment 

The vacuum alignment problem of the model can be solved by the supersymmetric driving 
fields method introduced in Ref. [33]. This approach exploits a continuous U{1)r symmetry 
under which matter fields have -R = +1, while Higgses and fiavon fields have R = 0. Such 
a symmetry will be eventually broken down to the R-parity by small SUSY breaking effects 
which can be neglected in the first approximation in our analysis. The spontaneous breaking 
of Ti3 can be employed by introducing the so-called driving fields with R = 2, which enter 
linearly into the superpotential. Five driving fields P^^ ^ind C,^ are introduced in 

our model, their transformation rules under the fiavor symmetry are shown in Table |2j We 
note that the driving fields and 6^ are necessary to stabilize the vacuum alignment under 
subleading corrections. At LO, the most general superpotential dependent on the driving 
fields, which is invariant under the fiavor symmetry group T13 x Z4 x Z2, is given by 

+ 92{v\^v)s.)u + hemu (27) 

In the SUSY limit, the vacuum configuration is determined by the vanishing of the derivative 
of with respect to each component of the driving fields 

dwy 



dx\ 
dwv 



fixl + /2X26 = (28a) 

fixl + /2X36 = (28b) 

fixl + f2Xi^3 = (28c) 

fsiXi^i + u\2i2 + ujxziz) = (29) 

f^iXiii + ^X2i2 + c^'xs^) = (30) 



8 



The above equations are satisfied by the ahgnment 





(X) = ^x| 1 h (0 = ^H 1 I (31) 
with the relation 

= —^v^, undetermined (32) 
Ji 

Without assuming any fine-tuning among the parameters /i and /2, the VEVs and 
are expected to be of the same order of magnitude, this is consistent with the conclusion 
drew from the charged lepton mass hierarchies. We note that if one component of x or ,^ has 
vanishing VEV, e.g. (.^i) = 0, Eqs. fl28al) - fl28cp imply (xi) = (X2) = (Xs) = 0. This means 
that the VEV of any component of the fiavons x ^ should be non-zero in order to obtain 
a non-trivial vacuum configuration. As has been shown in the previous section, at LO the 
Ti3 fiavor symmetry is spontaneously broken by the VEVs of and 77 in the neutrino sector, 
their vacuum configurations are determined by 



dwy 
dr]2 

dWy 

dw„ 



2fi'i?72?73 + = (33a) 

2fi'i?7i?73 + 5'20i?72 = (33b) 

2fi'i?7i?72 + fi'202?73 = (33c) 

Mei0i + 602 + e303) = O (34) 



The first three equations Eq.f l33ap -f l33cl) lead to two un-equivalent vacuum configurations §. 
the first is 





{<P)=v^\ 1 I , iv) =v,\ 1 I (37) 
with 

Vr, = —-^Vfj,, undetermined (38) 



29 



The second solution is 

(0) =\ \, {V)=\ V, I (39) 





Another one is 



^We note that the equations can be satisfied by two additional solutions as well. One is 

v<p^\ / 

(0) - I , (^) - I (35) 

/ V ^'»? 

(0) ^\ \, (^) = ( I (36) 

where w^, , v^^ and are undetermined. However, the above two solutions can be obtained by acting on 
the vacuum Eg.lfM]) with the elements T and respectively. Therefore these two solutions are equivalent 
to the configuration in Ea. (p9)) . 
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where f^, v^^ and v^^ are constrained. Using the ahgnment of x iii Eq.( pT]) . for the first 
solution shown in Eq.f l37p . we can immediately infer from Eq.([3i 



v^v^ = (40) 

We are led to the trivial solutions = = or t",^ = = 0, which can be removed by the 
interplay of radiative corrections to the scalar potential and soft SUSY breaking terms for the 
flavon fields. Therefore we choose the second solution in this work, this vacuum configuration 
can produce the results in the previous section. Then the minimization equation Eq. (!34|l 
implies 

^<^2 + ^<^3 = (41) 

This indicates that (02) and (^3) have to be equal up to a relative sign, thus is fully aligned 
as 



{^) = vA 1 (42) 




Starting from the vacuum configurations given in Eq. fjST]) . Eq.f p9|) and Eq. fl42l) and acting 
on them with the elements of the flavor symmetry group T13, we can generate other minima 
of the scalar potential. However, these new minima are physically equivalent to the original 
one, it is not restrictive to analyze the model by choosing the vacuum in Eqs. fl31|39lH2]) as 
local minimum. It is important to check the stability of the LO vacuum configuration, if we 
introduce small perturbations to the VEVs of the fiavon fields as follows. 





(0) = ^0 1 h iv) = v,\ 1 (43) 



After some straightforward algebra, we find that the only solution to the minimization 
equations is 

Xi = X2 = X3 = 0, y2= 1/3 = 

zi = Z3 = 0, wi=W3 = (44) 

Therefore the LO vacuum alignment is stable, then we turn to consider the magnitudes of 
fiavon VEVs. Since the VEVs of and 77 are closely related with each other through the 
equations Eqs. fl33ap -f l33cp . and they have the same charges under the auxiliary symmetry 
Z4 X Z2, we expect a common order of magnitude for the VEVs and f^. However, the 
VEVs of $^ = {XjO and $1/ = {(f), r]} can be in principle different and they are subject to 
phenomenological constraints. As we have shown in section 3.1, ($£) is responsible for the 
charged lepton mass hierarchies, and it is required to satisfy 

^ = ^ ~ ? - (45) 
A A ^ 

Among the three neutrino mixing angles, the solar neutrino mixing angle 6*12 is measured 
most precisely so far, the experimentally allowed departures of 6*12 from its TB value sin^ 9i2 = 
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1/3 are at most of order [IH3]- It is well-known that the superpotentials Wi, w^^, wf,^^ 
and are corrected by higher dimensional operators in the expansion (please see section 5 
and Appendix B for detail), which mostly can be constructed by including the combination 
on top of each LO term, thus all the three mixing angles receive corrections of order 
($,^)^/A^ (please see section 5 for detail). Requiring that the mixing angles particular 612 
lie in the ranges allowed by neutrino oscillation data, we obtain the condition 

^'-f~^<A. (46) 

The same condition follows from the requirement that the generated charged lepton mass 
hierarchies should be stable under subleading corrections. As a result, we can tolerate 
a moderate hierarchy between e and e' because of the strong constraint of the auxiliary 
symmetry Z^x Z2. It is a general conclusion that a hierarchy between the VEVs of the flavon 
fields can be accommodated in a "fully" separated scalar potential. This type of vacuum 
alignment is usually constructed to generate a large reactor angle [MIESI, i-e. 6*13 ~ Ac, 
although it is predicted to be exactly zero at LO. However, the subleading corrections to ^13 
turn out to be of order A^ in our model, as we shall demonstrate in next section. 



5 Subleading corrections 

It is crucial to guarantee that the successful LO predictions are not spoiled by subleading 
corrections, we will discuss this important issue in detail. The superpotentials we, w^^ , wl^^ 
and are corrected by higher dimensional operators, which arise from adding the products 
^u^ui invariant combination under Z4 x Z2, on top of the LO terms. Then the residual 
Z3 symmetry in the charged lepton sector would be broken completely by the subleading 
contributions. The lepton masses and mixing matrices are corrected by both the shift of the 
vacuum configuration and the higher dimensional operators in the Yukawa superpotentials. 
As a result, the mass matrices with subleading corrections can be obtained by inserting the 
modified vacuum alignment into the LO Yukawa operators plus the contributions of the new 
higher dimensional operators evaluated with the unperturbed VEVs. 

The subleading corrections to the vacuum alignment are discussed in Appendix B in 
detail. The inclusion of the higher dimensional operators in the driving superpotential 
results in a shift of the VEVs of the flavon fields, the vacuum configuration is modified into 

(V) = ^'^ (47) 

where v^, and remain undetermined, and all the shifts are of order e'"^ with respect to 
the LO VEVs, as is shown in Appendix B. Moreover, all components of (x), (0^ (0) iv) 
receive different corrections so that the LO alignment is tilted. 



(0) 











^x 














-Vd> 
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5.1 Corrections to the charged lepton mass matrix 

In the charged lepton sector, is corrected by the following operators 

e'l^l^lhd/K\ /i'=£$2$2/i^/A4, rm$2/,^/A3 (43) 

where all possible contractions among fields are understood. After lengthy and tedious 
calculations, we find that each element of charged lepton mass matrix gets a small correction. 
Concretely the corrections to the e row, /i row and r row are of order e^e'^f^, e^e'^w^ and 
ee'^frf respectively. As a result, the charged lepton mass matrix with subleading corrections 
can be parameterized as 

/ Ves'^ Ves' VeS^ \ / a{^e' al^e' ai^e^ \ 

me = y^e u'^y^e ujy^e \evd+\ a^^e uj^a^^e uja^^e ee''^Vd (49) 

V yr UJVt Uj'^Vt J \ ^0.32 ^^«33 / 

where the first term denotes the LO contributions, and the second term represents the 
corrections induced by the higher dimensional operators in Eq. fHHj) . The coefficients a^j{i,j = 
1,2,3) are complex numbers with absolute value of order one, their specific values are not 
determined by the fiavor symmetry. Furthermore, we have to consider the corrections from 
the shifted vacuum alignment. Since the shifts 6v-^. and 6v^. are of order e"v^ and e"v^ 
respectively, and the corrections to each matrix element contain one additional power of 
6v^./v^ or 6v^Jv^. Consequently, including these corrections only amounts to a redefinition 
of the afj parameter in Eq. (l49|) . As a result, the unitary matrix Ue%, which corresponds to 
the transformation of the charged leptons used to diagonalized me, is modified into 



Ue = ^\ 1 u cu^ \u'e (50) 




where U'^ is given by 



1 {Aee'^y 

U', = I -Aes'^ 1 {Ces'^r I (51) 



'Bes''^ —Ces''^ 




with 



Ae = (021 + w^a22 + wa23)/(37//,) 
Be = (4i + ^^2 + ^^aD/i^Vr) 

Ce = (4i + w'42 + uai^)/{3yr) (52) 

The charged lepton masses are corrected by terms of relative order e'^ with respect to LO 
result, therefore the charged lepton mass hierarchies predicted at LO are not spoiled by 
subleading corrections. 



'^Ue is the unitary matrix which diagonahzcs the hermitian matrix mime. 
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5.2 Corrections to the neutrino mass matrix 



First of all we focus on the corrections to the right-handed Majorana neutrino mass. We 
note that the modified vacuum alignment doesn't affect the Majorana mass at all, since flavon 
fields are not involved in the LO Majorana mass terms. The subleading corrections from 
higher dimensional operators are of the form ufUj^^/A^, thus every entry of right-handed 
Majorana neutrino mass matrix receives corrections of order e'^A instead of e'^A, which can 
be safely neglected. Then we move to consider the corrections to the neutrino Dirac mass. 
Among the independent terms of the type i^f£$|^/iu/A^, only the operators {h'fi(l)^hu)ii/A^ 
give non-zero contributions. As a consequence, the first and the third columns of Dirac 
mass matrix receive corrections of order e'^Vu- In addition to this correction, inserting the 
VEV shifts in the LO operators introduces independent corrections of order e'^Vu to the first 
and second column elements of the Dirac mass matrix. Including the above two kinds of 
corrections mentioned, we conclude that all the elements of Dirac mass matrix are corrected 
by terms of 0{e'^Vu)- With these results, we find that each entry of m'^^ except the (11) 
element receives corrections of order 5'^w^/A. Now we discuss the corrections to the Weinberg 
operators. The superpotential wf,-^-^ is corrected by the contraction 

{ihJK<l>t)iJA'' (53) 

Taking into account the contributions of the modified vacuum alignment in addition, we 
find all the elements of m^-^'^ receive corrections of order e'^Vu/ A. As a result, the overall 
correction to the light neutrino mass matrix is a most general symmetric matrix of order 
e"^Vul A. The neutrino mass matrix including subleading corrections can be written as 

'^11 '^12 '^13 \ 2 

m, = 5'2| V Z + «12 «22 «23 X ^^^^ 

'^23 '^33 / 

where the parameters x, y and z can be easily reconstructed from the LO couplings in Eq.( lT6l) . 
and the coefficients a^^ are OiX) unspecified constants. The matrix vny is diagonalized by 
the unitary transformation 

/ 1 \ 



J_ J_ 



K (55) 



where U'^ is close to an identity matrix with small corrections on off-diagonal elements, it is 
given by 



with 



\J[ = -{A,e"y 1 C,e'' (56) 





A, 



(r + ^-)(ar2 + <3)+^K2 + «r3) 

\/2[|x|2 - \y + z\^] 
z(r + z*){a^,, - + i{y - z){a^^^ - g-) 
—4(2/2* + 



^ _ - Q13) + - ^)(Qi2 - Q^3; .57^) 

^\\y-z\^-\x\^\ 
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The PMNS matrix is U, 
modified as 



PMNS — 



UjUiy, then the parameters of the lepton mixing matrix are 



sin^ e 



sin^ 9- 



sin 6'i3 = — 



23 



12 




(58) 



We see that all the three mixing angles receive corrections of order e'^ from both the neutrino 
and the charged lepton sectors. As is pointed out in section 4, the data on solar neutrino 
mixing angle 6*12 constrain e' < Ac- Then the reactor angle 6*13 is of order A^, it is within the 
sensitivity of the experiments which are now in preparation and will take data in the near 
future [3S1EZ]- Since three complex parameters which are related with three light neutrino 
masses are involved at LO, the light neutrino mass spectrum can be normal hierarchy or 
inverted hierarchy, and the phenomenological predictions of the model are just the generic 
results of neutrino mass matrix with TB mixing, e.g., degenerate neutrino mass spectrum is 
disfavored since strong fine-tuning is required to produce the observed mass squared differ- 
ences Am^^i and Am^^^, and the Oi/2/9-decay mass \mee\ in inverted hierarchy is generally 
larger than that in normal hierarchy. 

6 Phenomenological implications 

In the following, we shall investigate the physical consequences of our model, and the 
corresponding predictions are presented. We perform a numerical analysis by treating all the 
LO and NLO coefficients as random complex numbers with absolute value between 1/3 and 
3, the expansion parameter e and e' are set to the indicative values 0.05 and 0.22 respectively. 
In FigUl we plot the effective Oz/2/3-decay mass |mee| as a function of the lightest neutrino 
mass. The constraints which have been imposed to draw the points are the experimental 
values at 3a for the neutrino oscillation parameters Am^^^, Am^^;, sin^6'i2, sin^ 6'23 and 
sin^ 6*13 [TH3]. We also show the future sensitivity on the lightest neutrino mass of 0.2 eV 
from the KATRIN experiment [38] , and the horizontal lines represent the sensitivities of the 
future Oz/2/3-decay experiments CUORE [39] and Majorana [10]/GERDA III [H], which are 
15 meV and 20 meV respectively. It is obvious that the effective mass |mee| of inverted 
hierarchy (IH) is generally larger than that of normal hierarchy (NH). Since the bulk of data 
are predicted to be above the sensitivity of CUORE experiment for IH, the rare process 
Oi/2/3-decay should be observable in future, if the neutrino spectrum is IH. We note that 
most of the points fall into the region where the lightest neutrino mass is smaller than 0.02 
eV for NH spectrum, and a large set of points lie in the region of the lightest neutrino mass 
between 0.01 eV and 0.04 eV for IH case. The values beyond these regions, in particular the 
region of degenerate spectrum, are strongly disfavored. 

Finally, we show the sum of light neutrino mass as a function of the lightest neutrino 
mass in Figj2l The vertical line denotes the future sensitivity of KATRIN experiment, and 
the horizontal lines are the cosmological bounds [12]. The first one is at 0.60 eV, which is 
obtained by combining the data in Ref. [13], and the second one at 0.19 eV corresponds to 
all the previous data combined to the small scale primordial spectrum from Lyman-alpha 
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Figure 1: The effective mass |mee| as a function of the lightest neutrino mass. The red corresponds to the 
inverted hierarchy neutrino mass spectrum, and the black corresponds to the normal hierarchy case. The 
future sensitivity of 0.2 eV of KATRIN experiment is shown by the vertical solid line, the future expected 
boimds on \mee\ from CUORE and Majorana/GERDA III experiments are represented by horizontal lines. 

(Lya) forest clouds [B]. We see that the current cosmological information on the sum of the 
neutrino masses can hardly distinguish the NH spectrum from the IH spectrum. However, 
such a discrimination could be possible if these bounds are improved in the near future. 

7 Conclusions and discussions 

In this work, we have presented a T13 model for TB mixing based on the flavor symmetry 
Ti3 X Z4 X Z2. Both the charged lepton singlets e'^, fi'^, t'^ and the right-handed neutrinos 
vf are assigned as T13 singlets in this work. The light neutrino masses are generated as 
a combination of type I see-saw mechanism and Weinberg operators, and neutrino mass 
spectrum can be normal hierarchy or inverted hierarchy. In the charged lepton sector, the 
flavon fields $^ = {x? O break the T13 group into the Z3 subgroup at LO, and the symmetry 
breaking parameter e = (x)/-^ ~ (0/-^ controls the charged lepton mass hierarchies without 
invoking a Froggatt-Nielsen f/(l) symmetry. In the neutrino sector, the T13 group is entirely 
broken by the fiavon fields = {4>,ri} at LO, the symmetry breaking parameter e' = 
{4>)/A ~ (^)/A can be chosen to be of the order of Cabibbo angle Ac without spoiling the 
LO predictions and vacuum alignment. It is a noticeable feature that we can still reproduce 
the TB mixing even if the flavor symmetry is broken to nothing in the neutrino sector at 
LO. 

The subleading corrections are discussed in detail. The subleading operators contributing 
to lepton mass and vacuum alignment are obtained by inserting $^ into the LO operators 



T 1 1 I I I I I I 1 1 1 I I I I I I 1 1 1 — I I I I I I 1 1 1 I I I I I I 



Majorana/GERDA III 
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Figure 2: The effective mass |mee| as a function of the lightest neutrino mass. The red corresponds to the 
inverted hierarchy neutrino mass spectrum, and the black corresponds to the normal hierarchy case. The 
vertical solid line represents the future sensitivity of 0.2 eV from the KATRIN experiment, and the horizontal 
lines refer to the cosmological bounds. 

in all possible ways and by extracting the T13 x Z^x Z2 invariants. We have showed that all 
the mixing angles receive corrections at the level of O(A^), in particular, the reactor angle 
6*13 is predicted to be within the reach of next generation neutrino oscillation experiments, 
although it is small. Furthermore, since the neutrino mass matrix depends on three unrelated 
complex parameters at LO, the phenomenological consequences of the model are the general 
results of neutrino mass matrix with TB mixing, there are no model-dependent peculiar 
predictions. 

In the end, we discuss whether we can extend the T13 flavor symmetry from the neutrino 
sector to the quark sector. The most naive way is to adopt for quarks the same classification 
scheme under T13 x Z4 x Z2 that we have used for leptons. With such an assignment, 
both up and down type quark mass matrices are diagonalized by the same unitary matrix 
Ue shown in Eq.f ll2p . as a consequence the CKM matrix is a unit matrix at LO, this is a 
good first order approximation. The off-diagonal elements of the CKM matrix arise when 
the subleading contributions are taken into account. As has been showed in section 5, the 
subleading corrections to the three quark mixing angles are of order A^, the resulting CKM 
matrix should have the same form of the unitary matrix U'^ given in Eq.(l5T]). Therefore it 
seems difficult to reproduce the quark mixing without introducing new ingredients in the 
symmetry breaking sector. Furthermore, there is another lack if we adopt for quark the 
same structure as that in the lepton sector, the resulting mass hierarchies among the up 
type quarks are not realistic, although it is a satisfactory result that the mass spectrums of 
down type quarks and charged leptons are predicted to have the same pattern. Since the 
top quark mass is of order of the electroweak symmetry breaking scale, it is much heavier 
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than the remaining quarks, it is natural to assign quarks as 2+1 representation instead of a 
triplet. In the context of U(2) flavor group, this assignment has been known to give realistic 
quark mixing matrix and mass hierarchies [15]. Inspired by this assignment, it is usually 
suggested to extend the flavor symmetry group, which can produce the neutrino TB mixing, 
to its double covering in order to give a coherent description of all fermion masses and 
flavor mixings. The flavor models based on T' and I' groups [I2lll6], which are the double 
covering groups of /I4 and A5 respectively, have been studied extensively. These models can 
really lead to a good description of the observed pattern of quark masses and mixing besides 
reproducing TB mixing (or the gold ratio mixing pattern) for neutrinos. Following the same 
logic, we expect the double covering group of T13 can simultaneously describe the lepton and 
quark sector very well. 
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Appendix A: Representation matrices and Clebsch-Gordan 
coefficients of T13 

The Ti3 group has 7 inequivalent irreducible representations li, I2, I3, 3i, 3i, 82 and 82. 
The representation matrices of the generators 5* and T in these representations are given in 
Eq.(jl]) and Eq.([5]). In the following, we present the representation matrices of all the group 
elements for the three dimensional representations. The explicit expressions of the elements 
in the 3i representation are 
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p^ 
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p9 
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p 
C7 : 5^ = I p3 

00 pV \00p 

For the remaining 3- dimensional representations 3i and 82, the matrices representing the 
elements of the group can be found from those just listed for the representations 3i and 82 
by performing complex conjugation. The above representation matrices can help us to see 
clearly how the T13 flavor symmetry is broken in model building. Starting from the above ex- 
plicit representation matrices, we can straightforwardly get the product decomposition rules 
of the Ti3 group. In the following we use to denote the elements of the first representation 
of the product and /3j to indicate those of the second representation. 

• I2 ® 3i = 3i 




3i ~ I J'a^2 I (59) 

I2 ® 3i = 3i 

3i ~ I a;2^/32 1 (60) 
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3i 



3i 



3i 



3i 



af3i 



a(3i 







"2/33 + "3/32 
3is ~ I as/Si + ai/33 

ai/32 + «2/3i 

«2/33 - «3/32 

ai/32 - a2/3i 




li © I2 © I3 © 32 © 32 

li ~ aiPi + ^2/32 + "3/33 

12 ~ ai/3i + wa2/32 + ^^"3/33 

13 ~ aiPi + w^q;2/32 + was/Ss 
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3i ® 32 = 3i © 32 © 32 



3i © 32 = 3i © 3i © 32 



3i © 3i = 3 



15 ty oiA 



3i © 32 = 3i © 3i © 32 



"2/31 

«3/32 I (73) 

«2/33 I (74) 
as/?! 



«3/33 

3i ~ I ai/3i I (75) 
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>1A 



"2/3; 



'2 



«3/32 

«i/33 I (76) 

"2/31 
«3/3l 

«i/32 I (77) 
"2/33 



«2/3l 

3i ~ I a3/32 I (7^ 



ai/3; 



3 



«i/3i 

3i ~ I a2/32 I (79) 



"3/3: 
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"2/33 

«3/3i I (80) 
ai/32 



«2/33 + «3/32 

a3/3i + ai/33 | (81) 
ai/32 + a;2/3i 

"2/33 - «3/32 

a3/3i-«i/33 I (82) 
ai/32 - a2/3i 

ai/3i 

«2/32 I (83) 
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ai/32 + tt2/3i 
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ai/32 - a2/3i 



3i ©3i 

li ~ a!i/3i + a2/32 + "3/33 

12 ~ tti/3i + wa2/32 + (^"^asPs 

13 ~ ai/3i + w^a2/32 + i^a3/33 

"2/33 
3i ~ I asPi 

ai/32 

"3/32 

3i ~ I aiP-s 



0^2/3 
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3i © 325 © 32A 

3i ~ I tts/^s I (9^ 



>2S 



*2A 




"2/33 + "3/32 

a3/3i + «i/33 I (99) 
ai/32 + «2/3i 

"2/33 - «3/32 

a3/3i-«i/33 I (100) 
ai/32 - oi2h 



Appendix B: Vacuum alignment beyond leading order 

In this appendix, we shall discuss the subleading corrections to the vacuum alignment 
induced by the higher dimensional operators. At the next level of approximation, the driving 
superpotential in Eq.( l27|) is modified into + Sw^. Due to the constraint imposed by 
Z4 X Z2 symmetry, the correction terms are suppressed by l/A^. Concretely is given by 

^34 1 11 1 1 28 1 19 

1=1 i=l j=l i=l i=l 

where q, r,, t,, e, and /cj are order one coefficients, {Of , Of ,Of ,Of , Of } denote the 
complete set of subleading contractions invariant under T13 x Z4 x Z2. 
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2)335 


of = x°((xx)3,,(w)32s)32, of = x'{{xxh,si<Pvh.h. 
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= X°((XX)3is(0^)3: 


J325 


of = xMxxh,si<l>vh.h., of = x'iixxYMKsh., 
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= x\{xx)3^{vv)3i] 


)327 


of = xMxxhMvh.h., of = /((eO3,,(00)3,j32, 


*-^10 


= X°((a)3,,(00)3. 


)32' 


of = X'mhjr^r^h.h., of = x'mhjHh.h., 


*-^13 


= X°((e03.s(</'^)3. 


)32) 
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of = /(xOi.(00)32, of = /(xOi.(0^)32, 
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= x°(xOi.(00)3., 


0X° 
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= X°((xO3.(0r/)3j 


'32 


(102) 



of = pMxx)s.si^vhJis, of = p'iixxhMsJu. of = p%{xx)-sMs.)is. 

of = p'mh.si<p<Ph.s)u, of = p'mh^sivvh.hs, of = p'mhMh.)is, 
of = p'mh.{Hh.)u, of = p'iixOMWis, of = p'mh.{HWis, 

oto = p'mh.ivvh.s)u, of = pMxO-3Ms.)u (103) 
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of = e\ms.,{<p<phJu, of = e\m-s.sivvh.)u, of = e'mhM<Ph.)u, 
of = 0'mh.{<l>vh.)i., of = ^°((xO3.(00)3ji., of = e'mh.{<l>vWi.. 
ofo = oMxOMs^sh., of, = e\{xihMvh.)u (104) 



of = r/°((00)3,.(00)3,,)3., of = r/°((00)3,(00)3j3.., of = r/°((00)3,,(0r/)3j3., 

of = v'm)sMvh.)s., of = v'mh.i<Pvh.h., of = v'mhMvh.)s.s, 

of = v'mhMvh.)s,., of = v'mh.sivv)sjs., of = v'm)s,sivvh.s)s., 
0f=v\ivvh.i<Pvh.)s., 0f=V%ivvhA<Pv)s.)s., 0f=V%ivv)s,si<Pvh.h., 

Oi3 = v\ivv)3iivv)si)32, oi^ = v°{{vv)siivv)32s)s2, oi^ = v\ixx)3isixx)3,s)s2, 

of = V^ixxh^ixxWs^s, of = V° iixx) Sis ixO 3^)3^, of = r/°((xx)32(x032)32S' 

of = VMXX)3MS.)3,A, of = V%iXX)3.sm3j3., of = v\iXX)3.sm3.)3., 

of = AixxhMhJs., of = AmsM03.)s., of = v'mhMO-s.h., 

of = v'ms.{xOu. of = v'm3.{x0i.. of = v'm3.{x0u, 

of = V'mh J^03.)3. (105) 



of = e{{X<P)3,si<Pv)-3.)u, of = e{iX<P)3,A<Pv)-3ju, of = ei{X<P)-3.iVV)3,s)u, 

of = e{{x<Ph.{<Pv)3.)i., of = enxvhMsju, of = e{{xv)3.ivv)3.s)u, 
of = em)3M)-3.)u. of = e{m3.i<Pv)-3.)u, of = emh^Ms^u, 
ofo = e{m-s.{H)s.)u, of = mvh.{vv)s.)u, of = ernhMs^sh. im 

The subleading contribution Sw^ modifies the LO VEVs, then the new vacuum configuration 
can be parameterized as 
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(107) 



where the shifts 5v^^, 5v^^ and 5vn2 have been absorbed into the undetermined parameters 
wg, and f^. Similar to section 4, the new vacua is obtained by searching for the zeros of 
the F-terms, i.e. the first derivative of + 5w^ with respect to the driving fields p^, 
6^, Iff and . By keeping only the terms linear in the shift 5v and neglecting the terms 
proportional to 5f/A, the minimization equations become 

2/ifxHi + /2^cH2 + /2^xHi + «i^x^?^jM^ = 

2/1 ^^xHa + /a^cHi + «3^^x^«^J/^^ = 

MM^'^xi + ^^H2 + ^Ha) + ^x(Hi + ^^H)] + (^AV^Virvl/A'^ = 
/4k(Hi + ^H2 + ^^Hs) + ^x(Hi + ^H)] + «5^x^«^5M^ = (108) 
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where the coefficients ai{i = 1 — 5) are hnear combinations of the subleading coefficients 

ai = 2civ^/v^ + csVy.Vr,/{v^v^) + (-4c9 + 2cio)v^/v^ + 2cuv^v^/{v^vl) + (-3c2i + cu)vr,/v^ 

02 = 2(C3 - C::,)v^V^r,/{v^V^) + 2(Cio - Cu)v^/V^ + Ci5V^v'^/{v^vl) + 3(C20 - C26) + C27v'^/vl 

as = 2(ci - C6)vjv^ + C7t;xf J/ (v^vl) + 2(ci2 - Ci3)t;5t;^/ (v^v^) + 3c2o - 2c32 + c^^v'^/vl 
a4 = 2a;rifxt;^/(t;5t;0) - (4r4 + re)!'^/!;^ + 2rr,v^v^^/ (v^vl) - rjv^Vr^/ (v^v^) - - rgv^j/v^ 
as = 2u\v^nJ i^^v^) - (4^4 + tQ)v^/v^ + 2tr,v^v^/ (v^vl) - tjv^v^/ (v^v^) - ts - tgVn/i^09) 

The equations Eg. (II 08 p are hnear in Sv^^{i = 1,2,3) and = 1,2), and can be solved 

straightforwardly by 
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^a, - 4a2 + ^a,)—f- - (11 + 19a;)a4T^^ + (8 + 19uj)a,—f- (110) 
t;^ I3/2A2 39/3A2 39/4A2 

From the above equations, we clearly see that all the shifts 5v^^/v-^, 5v-^^/v^, 6v^^/v^, 6v^-^/v^ 
and 5v^^/v^ are of order e'^. The minimization equations for 5v^-^, ^f^j, ^^771 and 5vr^^ are 

2giv^5v^^ - g2V^Svn^ + biv^v^^/A'^ = 

g2V^Sv^^ + b2Vr,vl/A^ = 

2giVnSvr,^ + g2V^Svn^ + h^iV^vl/ = 

h{v^5v^^ + t'^^f^a + v^5v^^) + hiV^vyk^ = (111) 
where the coefficients 6j(z = 1 — 4) are given by 

d = [2(2ei5 + ei6)v^ + 2{en + eisyv^v^ + 2(2e2o + 621 + 622)^^5^'! + (2e23 + 2e24 + 3e25)^^x^| 

+ 2e2svl]/{vr,vl) 
hi = (4ei + 2e2)v^/v,^ - 2e^v^/v^ + Ci^vyvl + d 
h = -ee - 67 + d 
bs = -263 - 66 + 67 + eiov'^/vl + d 

64 = (ki + k2)v^v^/{v^v^) (112) 
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The solutions to Eg. fill II) are given by 



Sv4>2 _ M _ _ 

_ {hi92V<i> - 2b39iVr,)v^ 

Vr, 4:9fv^^ + 9ivl A2 

SVria _ i'2bi9lVri + b392V<p)v^ 

Vr, A9lvl + 9lvl A2 



(113) 



Obviously ^v^^jv^, 5v^^/v^, 5tif^-^/vn and ^f^g/wj, are of order e'^ as well. As is shown in 
Eq. fllOSp . the subleading terms proportional to rf are of the structures 77°$^ or 77°$^, the 
contributions of the latter operator to the vacuum alignment are parameterized in terms of 
the parameter d in Eq.f lll2p . If we have a large VEV of $1, with ($,^)/A ~ Ac, then the 
structure 77°$^ is dominant. On the other hand, if the VEVs of and $^ are of the same 
order of magnitude, the contributions of the two type of operators are comparable. 
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